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ABSTRACT 

We demonstrate a method for quantitatively comparing 7-ray pulsar light curves with magnetosphere beam- 
ing models. With the Fermi LAT providing many pulsar discoveries and high quality pulsar light curves for 
the brighter objects, such comparison allows greatly improved constraints on the emission zone geometry and 
the magnetospheric physics. Here we apply the method to Fermi LAT light curves of a set of bright pulsars 
known since EGRET or before. We test three approximate models for the magnetosphere structure and two 
popular schemes for the location of the emission zone, the Two Pole Caustic (TPC) model and the Outer Gap 
(OG) model. We find that OG models and relatively physical B fields approximating force-free dipole magne- 
tospheres are preferred at high statistical significance. An application to the full LAT pulsar sample will allow 
us to follow the emission zone's evolution with pulsar spindown. 

Subject headings: gamma rays: stars - pulsars: general 
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1. INTRODUCTION 

With the successful launch of Fermi, formerly GLAST, the 
number and quality of 7-ray light curves available for com- 
parison with magnetosp heric models has substantially im- 
proved (lAbdo et alll2010|) . The quality and uniformity of the 
data make this a good time to down- select models for pulsed 
magnetospheric emission, by comparison with the observed 
curves. In the present paper we describe a method of attack 
for this process and apply it to the new improved light curves 
provided by Fermi for several well-known 7-ray pulsars. 

The most popular pulsar models postulate regions in the 
magnetosphere where the force-free perfect MHD (E-B = 0) 
conditions break down, and in response to the uncanceled 
rotation-induced EMF E « — r x Q x B/c, the charges 
re-arrange so that non-negligible fields penetrate these lower 
density 'gap' zones, causing particle acceleration and radia- 
tion. Traditionally, one approximates the magnetosphere by a 
vacuum dipole field. 

There are three locations commonly discussed for these 
magnetospheric gaps. The earliest pulsar models focused on 
acceleration at the foot points of the open field line zone, 
the so-called "polar caps". This emission should arise from 
within about a stellar radius (Daugherty & Harding 1996), 
and should therefore suffer strong attenuation from one pho- 
ton (7-B) pair creation. The corresponding hyper-exponential 
cutoffs are not seen in the LAT pulsar spectra, and so it has 
been concluded that the b ulk of the pulsar emission must 
come from higher altitudes (lAbdo et alll2010|) . The first high 
altitude model posits "Holloway" gaps above the null charge 
surface, extending toward the light cylinder at tlc = c/ft 
(the "outer gap" model, Cheng, Ho & Ruderman 1986; Ro- 
mani 1996). More recently, it has been argued that acceler- 
ation at the rims of the polar caps may extend to very high 
altitude (Muslimov and Harding 2004), effectively pushing 
the polar cap activity outward. In this 'slot gap '-type picture 
emission can occur over a large fraction of the boundary of the 
open zone. The geometrical realization of this model makes 
low altitude emission visible from one hemisphere and higher 
altitude (r > vnc) emission visible from the other. This is 
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the "two pole caustic" (TPC) model (Dyks and Rudak 2003) 
which is thus intermediate between the polar cap (PC) and 
outer gap (OG) pictures. 

In Watters et al. (2009), we created a library of light curves 
for these three emission locations in a dipole field and classi- 
fied the phases of the dominant peaks. The results were pre- 
sented as an 'atlas' of pulse properties which could be used to 
predict the pulse multiplicity and phase separation of a given 
model for a particular viewing geometry. While this is use- 
ful for an overview of the possible light curves and provides a 
convenient reference to read off the allowed angles for a given 
model, the complexity of the modeled beams meant that an 
appropriate phase separation could be found in several mod- 
els for any one pulsar. 

These computations implicitly assume static co-rotating 
charges only. Of course, the presence of radiating pair 
plasma in the magnetosphere requires some particle pro- 
duction, and so the vacuum plus co-rotating charge model 
will be approximate at best. Therefore an alternative at- 
tack is based on numerical solutions for magnetospheric 
charges, currents and fields, assuming that pair produc- 
tion is so robust that the force-free condition holds ev- 
erywhere dCon topolous. Kazana s~& Fendtl 119991: ISpitkovskyl 
2006; Bai & Spitkovskvl l2009bl) . Suchfilled magnetospheres 
lack the acceleration fields required to produce powerful 7- 
ray beams, so the truth must lie in between, with some regions 
experiencing charge starvation and departing from force-free. 

Here we wish to show how the actual light curves from the 
computations can be compared with the data. We demon- 
strate that the present data quality are sufficient to start com- 
paring predictions between different assumed magnetic field 
structures and emission zone locations. In particular, while in 
many cases several models can produce plausible matches to 
the light curves for some angle, when we restrict our attention 
to the geometries demanded by external constraints, one can 
often strongly exclude certain models. 

For our model set, we start by amending certain in- 
adequacies in the vac uum computations pointed out by 
iBai & Spitkovskyl (l2009a|). This ends up maki ng modest dif- 
ferences to the results in IWatters etal (|2009). We also in- 
troduce a simple model to illustrate the field structure pertur- 
bations resulting from open-zone currents. Comparing with 
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the Fermi full band light curves, we find that lower altitude 
TPC-type models have great difficulty producing acceptable 
matches for many pulsars for all of these field geometries (at 
least for the assumptions made here) and that OG-type models 
are statistically strongly preferred. Further, in fermz-quality 
data we find that the light curve matches can often be signifi- 
cantly improved by the magneto spheric current perturbations. 
In a following communication, we will apply such analysis 
to the full set of Fermi LAT pulsars, marry the light curve 
fitting code to a population synthesis model, and extract addi- 
tional information on the pulsar beam structure and its evolu- 
tion with age as the pulsar spins down. 

2. MODELING METHOD 

Our basic assumptions are that i) the stable, time-averaged 
pulse profile is caused by emission zones locked to a set of 
'open' field lines, ii) these field structures are dominated by 
the dipole component of the field anchored in the star, since 
for young pulsars this emission occurs at many R*, Hi) that 
a co-rotation charge density distributed throughout the bulk 
of the magnetosphere causes the plasma to follow the field 
lines and rotate with the star and iv) in the charge- starved 
gap zones the uncanceled field rapidly accelerates e + je— 
to highly relativistic energies so that these 7-ray producing 
particles also follow the magnetic field lines. 

As has been well established, starting with Morini (1983) 
and discussed explicitly by Romani & Yadigaroglu (1995), 
the combination of field line sweep back, co-rotation of the 
radiating particles and travel time across the light cylinder 
tends to pile up the emission into 'caus tics' in pulsar phase. 
As noted by iBai & Spitkovskvl (l2009a|) this reflects the sin- 
gularity of the Jacobian matrix mapping from the emission 
surface to the antenna pattern on the sky. The singularity nat- 
urally produces sharp pulses with emission from a range of 
magnetospheric altitudes, a property shared by all competi- 
tive magnetospheric emission models. The shar pness of the 
peaks of the pulsed emission (lAbdo et alll2009b|) implies that 
the 'width' of the radiating zone is modest, at least in a time 
averaged sense. However, the resolved rise to peaks of the 
observed 7-ray pulses suggest either a wide range of emis- 
sion height, a fundamental radiation pattern with wings, or a 
dispersion in the direction of the emission zone, again in a 
time-averaged sense. 

Once one chooses a prescription for the magnetic field 
structure, one can follow radiation from the active emission 
zone to a (</>, () observer plane, where (j) is the pulsar spin 
phase, with phase zero when the surface magnetic dipole axis 
passes closest to the observer line of sight. The observer line 
of sight is at viewing angle ( to the spin axis. This two di- 
mensional intensit y map is the ' antenna pattern' for the pul- 
sar radiation (IWatters et all 12009). Once one specifies the in- 
clination a of the magnetic dipole axis to the rotation axis 
one has a complete geometrical description. Of course, the 
choice of the radiating zone and the (energy-dependent) in- 
tensity weight through that e mission zone depend on the radi- 
ation physics (Ro mani|[l996|) . In practice, we use simplified 
assumptions that capture the basic behavior, and look for de- 
partures from the resulting predictions as guides to improved 
modeling of the physics. For example, it is o ften assumed 
that the gaps self-limit to a fixed potential drop ( Arons 2006) 
which means that the 7- ray lumin osity will be proportional 
to the open zone current (lHardinglll98l|) . One natural way to 
do this is to imagine radiation reaction-limited particles pass- 
ing through a charge- starved gap spanning a fraction w of the 



open zone of the pulsar magnetosphere. In this paper we take 
the emission surface to be relatively thin, lying a fraction w of 
the way from the last closed field line surface to the magnetic 
axis. For w not too close to one, this approximates the con- 
stant voltage limit if one takes the heuristic 7-ray luminosity 
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This implies an efficiency r\ = L 1 /Esd equal to the frac- 
tional gap width w. Realistic models must, of course, have a 
finite thickness for the emission zone. We find that this pro- 
duces relatively small perturbations to the light curve shape if 

Sw < w. 

The phase-averaged 7-ray flux for the Earth line-of-sight 
F }) S is related to the total 7-ray luminosity by 

L 7 = ^7rf Q F obs D 2 (2) 

where D is the distance to the pulsar ( Wa tters et alll2009l) and 
fn is the 'flux correction factor'; an isotropic emitter corre- 
sponds to = 1. Notice that we need this quantity to infer 
the actual pulsar 7-ray efficiency 
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for comparison with model p redictions . A presc ription for 
fn is given in Equation (4) of WattersetaJ (|2009); the solid 
angle corrections are applied correctly in the values in that 
paper (albeit for the 'Atlas' -style beam computation). 

2. 1 . Magnetic Field Structures 

To start, one must choose a field line structure. In the ear- 
liest l ight curve modeling (Morini 1983; IChiang & Romanil 
1 19941) static dipole magnetic field lines were assumed. Emis- 
sion was assumed beamed along these field lines, which were 
taken to occupy the rotating frame, and the resulting aberra- 
tions and time delay formed the observed pulse. While this 
"static" (Stat) model had the right topology (closed zone plus 
flaring open zone), and already predicted (for outer magneto- 
sphere emission) the generic double 7-ray pulse lagged from 
the low altitude radio pulsations, it clearly is not a physically 
consistent picture. 

The next step, introduced by iRomani and Ya digaroglu 
(1995) and used widely by other modelers, is to use the re- 
tarded dipole instead of the static case. This is commonly re- 
ferred to as the Deutsch (1955) field approximation. However, 
at large radius from the star, this is effectively just a rotat- 
ing point dipole which is most simply expressed by Kaburaki 
(1980) and can be re- written as 
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m(t r ) + rm(t r ) + r 2 fh(t r ) jr' 



+r (r • 3m(t r ) + 3rfn(t r ) + r 2 rh(t r ) /r 3 ) 

m{t r ) = m (since cosoot r x + sina sincj£ r y + cosa z) (4) 

with m the dipole magnetic moment evaluated at the re- 
tar ded time t r = t — r/c. This is the field structure used 
in (IWatters et al 2009) and so we refer to it as the "atlas" 
(Atl) field. It has the virtue of having the intuitively expected 
'sweepback' as one approaches the light cylinder (Figure [T]). 

This field is stationary in the non-inertial frame rotating 
with respect to the lab f rame, where particle s are assumed to 
follow the field lines. In IWatters et all (120091) and earlier com- 
putations by a number of groups, the particle velocity in the 
lab frame was incorrectly computed as a co-rotation induced 
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boost to a relativis t ic parti cle following the rotating field lines. 
iBai & Spitkovskyl (l2009a|) correctly point out that the particle 
velocities in the two frames are instead connected by a sim- 
ple coordinate transformation. These authors assumed that 
the earlier computations (incorrec tly) used the field (4 ) in the 
rotating frame. The boost error of IWatters et all (l2009|) is con- 
ceptually different, but mathematically equivalent to this as- 
sumption. In practice, the coordinate transformation allows 
one to compute the velocity /3y along the magnetic field line 
for the particle forced to co-rotate with the star at a highly 
relativistic lab frame velocity \/3q\ — > 1 

j3o=P\\B + (lx?/c (5) 

which gives 

/3f| = -B • (fi x f/c) + { [B ■ (fi x f/c)] 2 - (fix f/c) 2 + 1} 1 / 2 

(6) 

so that substituting into Equation (5) gives the direction of the 
particle as viewed in the lab frame. Note that /3n can become 

very small for particles traveling nearly along fi x f for r± ~ 
tlc, i-e. nearly tangent to the light cylinder in the 'forward' 
direction. In this case small effects of particle inertia or field 
line perturbations will 'break open' such field lines and the 
particles should leave the light cylinder and contribute little to 
the pulse emission. Although we do not follow these physical 
effects, such field lines do not affect our light curves, since we 
taper the emissivity (see §2.4 below). 

In any case, to have the plasma static in the rotating frame 
requires the presence of the co-rotation charge density that 
produces the lab frame field 

E = (fi x f) x B/c. (7) 

This is assumed present in the 'Atlas' and similar computa- 
tions, and is similarly assumed here. We refer to the revised 
computation as the 'pseudo force-free' (PFF) case, since it 
contains co-rotation enforcing charges, but no currents. These 
PFF computations, for both TPC and OG-type models, give 
light curves slightly distorted from in the 'Atlas '-type com- 
putations and appear to m atch well to the sample an tenna 
patterns and light curves of Bai & Spitko vskyl (1200 9a). The 
principal effect, as noted by these authors, is that for small 
w < 0.02 the low altitude emission that dominates the second 
pulse of the TPC model has decreased phase folding, making 
the Jacobian non- singular. There is a pulse peak but with- 
out a true caustic there is no sharp pulse. This also weakens 
somewhat the first TPC peak for such high E pulsars. For 
pulsars with larger w and for the higher altitude emission of 
the OG model there are departures from the Atlas '-style com- 
putation, but the effects are more subtle. In the Appendix we 
present arrays of example light curves for the TPC and OG 
models. These may be compared with, and considered as up- 
dates to, the similar light curve panels in Watters et al (2009). 
We also update the figures describing the pulse width and flux 
correction factors and their variation with viewing geometry, 
for comparison with that paper. 

2.2. Current-induced B Perturbations 

The gap-accelerated radiating charges alone produce some 
currents in the open zone, and we expect the gap-closing pair 
front to produce densities comparable to the co-rotation value 
for at least some of the open field lines. This current will 
perturb the magnetic field of Equations (4). 
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Fig. 1 . — Illustrations of the basic dipole field geometry. Left: Static, Right: 
retarded point dipole (as used in the 'AtT and 'PFF' computations). Top - 
the position of the open zone boundary (polar cap boundary) in magnetic 
coordinates 0b(4>b) for magnetic inclination a = 70°. Middle - magnetic 
field lines in the equatorial plane for a = 90°. Bottom - field lines for 
a = 80° above one polar cap with lighter (blue) showing the lines which 
cross rjsrc, the 'Outer Gap' lines. 



In the spirit of attempting analytic amendments to the vac- 
uum model, we applied the approx i mate perturbation field 
computed by Musl imov and Harding ( 2009) for a pair-starved 
current flow in the open zone. This is expressed in magnetic 
coordinates (r#, #b, (/>#), where the perturbation amplitude 
e' determines B' — e' [2m fi/(r 2 c)] with m the magnitude of 
the star's dipole moment 

and s = cos rot = cos a cos 0b + sin a sin b cos <p b . This 
perturbation field is defined with respect to the static vac- 
uum dipole. We take e' positive for outward- flowing electrons 
(i.e. negative current) above the magnetic pole passing clos- 
est to the Earth line-of- sight. This pole is generally inferred 
to produce the low-altitude radio pulsar emission. If this pole 
has the opposite current (i.e. opposite magnetic polarity) we 
have e' < 0. Of course an electrodynamically self-consistent 
model will have a particular sign of e' for, e.g. the positive 
magnetic pole, but our PFF model includes no current and so 
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is insensitive to the sign of B\ only this perturbation current 
breaks the degeneracy. 

To approximate the retarded solution, we matched this per- 
turbation to our magnetosphere, mapping B(vb, 0, 0#) to 
the magnetic axis of the swept-back dipole. The perturba- 
tion field becomes singular along the rotation axis where 
s — > ±1. By exponentially tapering this singularity (ex 
I _ e O|-i)/cr s ^ w ^ m _ o 25) we were able to integrate the 
field lines to determine the last closed field line surface (see 
below). Unfortunately, for lines near = 0, tt the resid- 
ual effects of this singularity still dominate, making the cap 
boundary at these azimuths unstable (Figure [2]). 

Accordingly we proceeded with a simpler toy perturbation 
field computed by integrating a constant line current passing 
through the star along the magnetic axis (the current extends 
to cylindrical radius r±_ = rs'mO = 1.2Rlc = 1.2c/ Q to 
minimize end effects). This perturbation field 

B'{r) (x I dlxr (9) 

JB(r b ,0,4> B ) 

is treated as static in the rotating frame and is summed with 
the field of (4). The effect of this simple perturbation field 
is dominated by the twisting of field lines along the magnetic 
axis, as can be seen by the shift of the 'notch' foot-points 
of the last closed field line ca p boundaries in Figure |2) This 
rotation is also present in the Muslimov and Harding ( 2009) 
field, but is dominated by the singularity effects except for 
nearly aligned rotators a « 0. For ease of comparison, we 
scale the perturbation amplitude for fields (8) and (9) to the 
underlying dipole field (4) so that 

e = (B'(0.5r LC , 0, 4>))/(\B(0.5r LC , 0, </>)\) 

where we average over a sphere at r = O.brLC- Again the 
sign of e is determined by the current; we expect opposite 
values for the two poles. 

None of these fields is fully realistic. In particular the 
current-induced B field perturbations used above are not com- 
plete consistent field/current systems. Nevertheless, given 
that the vacuum field geometry has already provided encour- 
aging successes in approximating the observed light curves, it 
is useful to explore how current systems can perturb the mod- 
eled beam patterns and pulse shapes. As an example, Figure 
[3] shows the 'antenna patterns' in the hemisphere containing 
emission from above the null charge surface (the OG case be- 
low). The unperturbed beam pattern (middle) suffers opposite 
distortions in the two hemispheres. 

2.3. Gap Locations 

For each field structure, we assume the appropriate co- 
rotation charge density so that the field pattern is stationary 
in the rotating frame. We define the open zone by identify- 
ing field lines tangent to the light cylinder (actually at slightly 
smaller radius for numerical expediency). These field lines 
are traced back to the surface where they define the open- 
/closed-zone boundary. In Figure [T] we show the foot-point 
locations, equatorial plane field lines (for a = 90°) and the 
last closed field line structure for an inclined a = 80° rota- 
tor. Our radiating surface is always in the open zone toward 
the magnetic axis from this boundary. This plot shows the 
static ('Stat') and point retarded dipole ('Atl' & 'PFF') field 
structures. The active field lines for the outer gap are lighter 
(plotted in blue). 

With the emitting field lines defined, we must choose 
the extent of the radiating surface. For the TPC model, 
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Fig. 2. — Open zone boundary (polar cap) foot-points: magnetic colat- 
titude b on the star surface as a function of magnetic azimuth 4>b for 
the a = 70° PFF fi e ld wi th current-induced perturbations. Top: the 
Muslimov and Harding ( 2009) pair-starved field, mapped to the retarded 
dipole. Bottom: retarded dipole with perturbations from magnetic axis 
line currents. Notice that the currents 'twist' the field lines for both mod- 
els, shifting the 'notch' at d> w 0.3. Although the singularity in the 
Muslimov and Harding ( 2009) field (upper panel) has been smoothed, we 
still see dramatic instability for lines directed near the rotation axis (4>b ~ 
0.4 -0.6 and 0s « 0.8 - 0.1). 




Fig. 3. — Antenna patterns for current-induced perturbations. Pulsar light 
curves are obtained from horizontal cuts across these images, after normaliza- 
tion by [sin£] _ 1 . Here we show emission from above the null charge surface 
(OG case) for a = 60°, w = 0.1. The middle panel shows the unperturbed 
PFF structure. The upper (e < 0) and lower (e > 0) panels show the effect of 
a perturbing line current along the magnetic axis. Note that the front caustic 
is better closed for e > 0, but for negative e the 'notch' structure is more pro- 
nounced. The dashed lines show the range where the viewing angle ( passes 
close to the low altitude emission from the opposite magnetic pole, i.e. when 
one expects a radio pulse. Measurement of light curve perturbations can thus 
determine the sign of B for a given magnetic pole. 

we follow the original definiti on (iDyks and Rudakl 120031: 
iDyks. Harding, and Rudak 2004), taking the emission surface 
to be the full set of last-closed field lines, but stopping emis- 
sion once the radial distance is r > vlc or once the dis- 
tance from the rotation axis exceeds r± > 0.75rz,c- This 
cut-off is required to avoid the higher altitude pulse compo- 
nents (as used in the OG picture), otherwise the modeled light 
c urves generally have three or four pulse components. As 
in lWatters et all (|2009|) we make this model more physical by 
placing the emission surface in the open zone, separated from 
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the last closed surface by the thickness w appropriate to the 
pulsar E (Equation 1). 

In the outer gap (OG) model radiation is produced on field 
lines crossing the null-charge surface • B(Rnc) = 0), 
with emission starting near this crossing point and extending 
toward the light cylinder. Again, we assume a gap thickness 
w. Field lines that do not cross the null charge surface before 
vlc do not radiate. 

The outer gap emission should be dominated by a section 
of the flared cone of field lines above a given pole. For these 
simple vacuum-based models, some field lines just inside the 
last closed surface extend for a very long distance before exit- 
ing the light cylinder or crossing the null charge surface. For 
example, some field lines arc over the star and cross B = 
for the first time far from the star with a path length (in units of 
Rlc) s ^> 1. These field lines define a disjoint, 'high altitude' 
second gap surface, that does not in most cases appear to be 
active (although the HF pulse components of the Crab might 
represent such emission; Moffet & Hankins 1996). Phys- 
ically, we expect that such lines arcing near Rlc will be 
opened in a real magnetosphere with currents and particle in- 
ertia. Also, gaps on field lines starting at large vnc will likely 
be inactive if the soft emission needed to close the gap arises 
fr om other field lines at much sm aller vnc- Accordingly, as 
in lRomani and Yadigaroglul (Il995|) . we apply a path-length s 
cutoff in the gap surface emissivity. Here the functional form 
used (for s > 1 + 2,s N c,min) is 

F(s) oc e -[(*- 2 *Nc, m i n -i)/a s ] 2 5 a s =0.1 (10) 

with 5NC,min the lowest null charge crossing for any active 
field line and all distances in units of Rlc- A detailed sur- 
face emissivity weighting awaits a physically realistic spectral 
emission model. 

While this work was being prepared for publication, a new 
discussion of light curve formation in fully force-free mod- 
els has been presented by Bai & Spitkovsky (2009b). These 
authors note the limitations with pure vacuum modeling men- 
tioned above. They also argue that in fully force-free mag- 
netospheres if the radiation is associated with the return cur- 
rent sheet, neither TPC nor OG pictures produce realistic light 
curves, and propose a new 'Annular Gap' (AG) geometry that 
has emission extending beyond the light cylinder. The foot- 
points of the AG field lines are the same as those used here for 
the TPC geometry (for a given w) but in iBai & Spitkovskyl 
(2009b) the co-rotating pulsed emission is taken to extend 
well beyond the the light cylinder. Thus, this model covers 
our OG zone plus additional emission at both high and low 
altitudes. As we see below, the vacuum dipole field (possibly 
with current-induced B perturbations) does a rather good job 
of matching the observed pulsations - it will be interesting 
to see if the force-free numerical fields using the AG or OG 
geometries can match the data comparably well. 

2.4. Computational Grid 

We identify the last closed field line surface (those field 
lines tangential to the speed of light cylinder) and trace these 
field lines to the polar cap surface at r* . We follow the field 
line structure to < 10 -4 rccs allowing good mapping even 
for pulsar periods of > Is. A single magnetosphere model 
can serve a range of pulsar periods, since the OG and TPC 
model features are computed as distances in fractions of vlc- 
The pulsar period only affects R*/Rlc and hence the inner 
boundary of the emission zone. 



Models are computed for all magnetic inclinations a. Pho- 
tons are emitted along the particle path in the lab frame (mov- 
ing tangent to the local field line in the co-rotating frame) and 
then phased, including light travel time across the magneto- 
sphere, to add to the skymap antenna pattern, as would be seen 
by distant observers. We assume uniform emissivity along all 
active field lines (except for the taper at large pathlength s as 
noted above). In practice we trace the trajectories for ~ 10 4 
photons per field line, which produces adequately smooth 
light curves for the present computation. Models are com- 
puted for a range of gap widths from w = 0.01 to w = 0.3, to 
accommodate a range of pulsar E. 

2.5. Lightcurve Fitting 

Our next task is to take this family of light curves and com- 
pare with the observed pulse profiles. We take our 7-ray data 
from the full ban d E > 0.1 GeV pulse profiles published in 
AbdoetaJ <2010). These use the first six months of Fermi 
data and give good statistics for the relatively bright pulsars 
considered here. In this Fermi pulsar catalog, a background 
level is estimated from the surrounding sky. Interestingly, 
some pulsars have significant steady (DC) emission. Match- 
ing this emission is an important test of the model - accord- 
ingly we fit to the total emission, pulsed a nd un-pulsed, afte r 
subtracting the background, as estimated in lAbdo etall (|2010). 
In particular, we do not fit the un-pulsed flux as a separate de- 
gree of freedom - the excess above the LAT pulsar catalog 
background is determined by the magnetosphere model. 

We take the measured pulsar E and compute w according 
to Equation (1), then use the closest value from the model 
grid. For this grid we take all a (in one degree intervals) 
and fit the model light curve for each a and £. To com- 
pare with the data we bin a given model curve into the phase 
bins (25 or 50/cyc le) defining the light curve for the pulsar in 
lAbdo et all (l2010|) . The total model counts are normalized to 
the background- subtracted data counts and we then search for 
the best-fit model to the combined (pulsar plus background 
emission) light curve. 
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Right) for Stat, Atl and PFF field structures. Top are OG models, bottom are 
TPC. The green lines show £ allowed by CXO PWN fitting. The ellipses give 
the combined constraints including radio polarization data. 

A variety of fitting statistics work well. Simple \ 2 fit- 
ting generally finds the best solutions, but is controlled by 
the more numerous off-pulse and 'bridge' emission phase 
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Fig. 5.— Pulse profiles for the TPC (left) and OG (right) models for Vela 
using the PFF field structure. The dotted pulse profiles show two periods of 
the LAT (E > 0. lGeV) Vela data from ?. The first model pulse in each panel 
is for the PWN/RVM-fit a, £, with no current perturbation while the second 
is for the parameters of the global OG fit minimum. The TPC fit quality 
varies only weakly with angle and does not improve substantially away from 
the RVM values. 

bins and so occasionally provides poor discrimination against 
models with unacceptably weak peaks. Choosing the model 
with the smallest maximum model-data difference for any bin 
Max [Abs (Modi — Di)] does an excellent job of matching the 
peaks, but less well on the bridge flux. We find that an expo- 
nentially tapered weighting of the largest model-data differ- 
ences is robust and produces sensible results. Thus we mini- 
mize 

_ (Modj ~ Di) -i(|Modi-Di|)/n 

Xn ~ Di 

with z(|Modi — Di |) the index of the phase bin differences, 
sorted large to small. Small n <C 1 just uses the biggest dif- 
ference bin (and thus focuses on the peaks), large n goes over 
to classical \ 2 weighting. The method is only weakly sensi- 
tive to n and an e-folding weighting n « 3 worked well. 

Our reference phase = occurs when the dipole axis, 
spin axis and Earth line of sight all in the same plane. It 
is generally assumed that radio emission arises at relatively 
low altitudes (but see Karastergiou and Johnston 2007), and 
that the radio pulse occupies a fraction of the open zone 
field lines at low altitudes. In lAbdo et al (|2010|) phase zero 
is taken to be the amplitude peak of the main radio pulse. 
However, the radio pulse profile often exhibits 'patchy' illu - 
mination of the emission zone dLyne and Manchester 1988). 
Thus the radio peak may not mark the phase of the closest 
approach of the magnetic axis. Additional geometrical in- 
formation comes from the linear polarization sweep of the 
radio pulse, with the maximum sweep rate associated with 
the line-of-sight passage past the magnetic axis. This max- 
imum often occurs well away from the intensity peak. Fi- 
nally, if the radio emission arises well above the star sur- 
face, relativistic effects can perturb the phases of the closest 
approach (ideally pulse intensity maximu m) and sweep rate. 
Blas kiewicz. Cordes & Wassermanl {1991) argue that the ap- 
proximate phase shifts are 

0int ~ ~r/r LC 0p o i ~ +3r/r LC . 

Thus polarization information can give additional constraints 
on the absolute phase, but altitude (plus mode changing and 
other sweep perturbations) can lead to substantial additional 
uncertainty. Accordingly, we have the option to fit light 
curves with the model phase free. Generally, we allow cf) 
to shift by as much as ±0.1 of a rotation, and retain the 
best fit model. When no radio pulse is known we allow full 



5(j) = ±0.5 freedom in the pulse phase. The phase for the 
best-fit model is recorded in the (a, Q map, as well as the 
model light curve and properties of the pulsed emission, for 
us e in further analysis. Recalling that the radio phase = 
in lAbdo etal (|2010|) is set at the pulse intensity peak, we ex- 
pect that our fit, which estimates 0o for the true magnetic axis, 
will tend to give small S(j) > for emission from finite pulse 
heights. However, we expect that 5<j) will be less than 0p o i. In 
practice, patchy pulses and polarization sweep uncertainties 
make S(j) of any modest amplitude plausible. 

In considering the values of the fit statistic % n , we must re- 
mind the reader that these are not complete models. Detailed 
radiation physics will certainly change the light curve shape 
from the uniform weighting assumed here. Also, the lack of 
self-consistent magneto spheric currents must, at some level, 
cause light curve shape errors. This, coupled with the excel- 
lent signal-to-noise of the Fermi light curves guarantees that 
these will not be statistically "good" fits. Monte Carlo exper- 
iments with the Xn fit statistic show that for Poisson realiza- 
tions of a 50 bin light curve, the typical value is (%3) = 7.5 
while 99% and 99.9% values are 15.5 and 19. Thus fit dif- 
ferences of A%3 ~ 8 are statistically significant at the nor- 
mal 2.5<j level. This gives a guide to interpreting the relative 
goodness of fits to the various models. 

3. FIT EXAMPLES 

We start with the archetype 7-ray pulsar, Vela. For this pul- 
sar we show the fit statistic \ surfaces for the three approxi- 
mate dipole field structures (Stat, Atl, PFF) discussed above, 
with darker colors indicating better model fits. The upper row 
is for OG models, the lower row for TPC. For OG all field 
structures show a similar topology, with best-fit solutions at 
C ~ 70 - 80° and a w 60 - 90°, and no solutions at small 
a, C The swept-back models (Atl and PFF) have a tail of so- 
lutions to smaller a. The TPC scenario can produce models 
for nearly all a, (. How ever, except for the 'Atlas' geometry 
which, as also noted by iBai & Spitkovskyl (l2009a|) . has un- 
realistically sharp caustic pulses for small gap width w, the 
models produce much poorer fits than the OG scenario, as 
they have weak pulses and far too much off-pulse (DC) emis- 
sion. 

3.1. External Angle Constraints 

The model fits in (a, Q space are especially useful when 
we have external constraints on these angles. Some care is 
needed, however, in applying angle constraints from other 
wavebands. For Vela and a number of other young pul- 
sars with bright X-ray pulsar wind nebulae (PWNe) we have 
been able to place relatively model-independent constraints 
on viewing angle by fitting for the i nclination C of the rela - 
tivistic Doppler-boosted PWN torus (iNg and Romanil 12008). 
This method does not, however, distinguish the sign of the 
spin axis. Therefore the X-ray torus fits allow two possi- 
ble viewing angles between the positive spin axis and the 
Earth line-of-sight C = C and (' = 180° - C In fact, in 
INg and Romanil (|2008) the torus position angle was restricted 
to be 0° < \£ < 180°, so the angles quoted in this paper may 
represent either ( or (' . 

In the context of the rotating vector model 
(Radhakrishna n and Cookd Il969|) radio polarization data 
can also constrain the viewing angles. In most cases, the 
small range of phase illuminated by the radio pulse allows 
only an estimate of the magnetic impact parameter 

(3 = (-a « sm- 1 [sina/(d^/d0) max ] (11) 
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where the maximum rate of the polarization position angle 
sweep \I/((/>) occurs at <j) po i — 0, near the closest approach 
to the magnetic axis. Here the sign of the sweep is mean- 
ingful, determining whether the line of sight is closer to or 
farther from the positive rotation axis than the observed mag- 
netic pole (at inclination a). Thus, when we have an inde- 
pendent PWN estimate of (, we have two possible a values 
- e.g. for a negative sweep we get f3 < and a = ( — f3 or 
a = (' — f3. In this paper, we tabulate the possible values of ( 
and but save space by plotting only in the positive rotation 
hemisphere. Thus there are two values possible for a from 
Equation (10), but when a or ( > 90°, we reflect back into 
the positive hemisphere for plotting purposes. 

Occasionally, when the radio pulse is very broad or when 
the pulse profile presents an interpulse, the radio polarization 
can make meaningful estimates of both a and (, from fits to 
the full polarization sweep 



tan(^ + ^ ) 



smft sin 



sm£ cosa — cos£ sina cos((/> — 0o) 

where the pol arization has the absolute position angle \I>o at 
i)Q. iKeith et all (|2009) have recently presented sev- 



<PPol 

eral examples. When the radio illumination is insufficient for 
a full solution, we can still break the degeneracy of the fits 
using the X-ray ( constraints. Finally, if the magnetic incli- 
nations orientations are, at least initially, isotropic we expect 
the inclination choice giving a closest to 90° to be preferred 
on statistical grounds. 

In Figure [4] t he tw o possible RVM a solutions for Vela 
(John ston, et all l2006|) consistent with the PWN data are 
shown by the green circles. Note that the OG and TPC mod- 
els have minima close to the solution at larger a (bold circle). 
In general, with constrained a, £ values one can make a dis- 
crimination between the models. Note that for the OG model 
the position of the minimum is closest to the preferred angle 
for the relatively physical 'PFF' case. The values of the fit 
statistic are substantially worse for TPC at all locations near 
these preferred angles (Table 1). Figure [5] shows the PFF light 
curves at the RVM angles, and at the global fit minimum, for 
the TPC and OG scenarios. 

3.2. Testing Field Perturbations 

The best-fit models in Figure |4] are several a off of the 
PWN/RVM angles. This suggests either that there are sys- 
tematic errors in these angle estimates or that the true mag- 
netosphere structure is slightly different than the PFF esti- 
mate above. The measurement of such offsets for a num- 
ber of pulsars should enable us to map the required geome- 



try differences and test physical models for their origin. As 
an example, here we illustrate the perturbing effect of mag- 
netospheric currents. Figure [6] shows the x plane for fits in- 
cluding a magnetic axis current (Equation 9). The current 
shifts the fit minimum and distorts the light curves. For a 
best fit in the RVM-determined region positive currents are 
required for the OG model. The best models within ~ la of 
the PWN+RVM angles also prefer a small positive current. 
For negative e the best fit region shifts further from the ex- 
ternally determined angles, although the match to the pulse 
shape is quite good, reaching a global minimum (x = 34) 
at (a, C) = (82°, 71°) for e = -0.2. Since the global min- 
imum does not shift precisely to the PWN/RVM angles, we 
infer that (unsurprisingly) the true field geometry is only ap- 
proximated by this simple model. The TPC model fits are not 
improved by including finite e; the best minimum anywhere 
near the external PWN/RVM angles has x ~ 200. 

To summarize, for Vela with the external constraint, PFF is 
clearly preferred and the best fit models are relatively close 
to the known a, (. The rotation axis is inferred to be ori- 
ented out of the plane of the sky, we are viewing at ( w 65° 
and the magnetic inclination is large at a w 75°. Positive- 
current field perturbations are preferred. The 7-ray pulse fits 
imply that 0o lies earlier than the radio peak (and earlier than 
the maximum sweep rate). This would suggest that Vela has 
patchy radio emission dominated by the 'trailing' side of the 
radio zone. 
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Fig. 7.— Goodness of fit surface for PFF models of PSR J1952+3252, 
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indicated by the maximum polarization sweep rate d^/d</>| m ax are shown 
as two green wedges. 
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3.3. Other EGRET Pulsars 

The next most energetic pulsar in our test set is PSR 
Jl 952+3253. As for Vela, the three field structures give quite 
similar %-plane maps. In (Figure U} we show the x maps for 
the PFF model with three values of B' . This pulsar has two 
narrow peaks separated by A = 0.49 in phase, so we find 
that the TPC model can produce reasonable solutions, as indi- 
cated by the relatively dark regions in the lower panels. This 
object is presently interacting with the shell of its supernova 
remnant, CTB80, and so produces a PWN bow shock rather 
than a torus, precluding a ( measurement. However from the 
max sweep rate of the radio polarization from Weisberg et al. 
(1999), we can infer a constraint on the magnetic impact pa- 
rameter (two green diagonal bands). The sweep maximum 
appears to lie towards the tail of the broad radio pulse, consis- 
tent with our fit S(j) w 0.05; again the radio pulse peak would 
represent patchy emission, in this case from the leading edge 
of the emission region at relatively low altitudes. 

Both models prefer relatively large a and (, and each can 
produce acceptable light curves within the region allowed by 
the radio data. If improved radio measurements or other an- 
gle constraints can be obtained, it may be possible to select 
between these two models for this object and also measure 
the altitude of the radio emission zone. 

Next and only slightly less energetic is PSR B 1706— 44 
= PSR J1709-4429. This object, in contrast to PSR 
J1952+3252, has a narrow double pulse with A = 0.25. 
Here OG models clearly prefer the small a RVM solution. 
While the best OG fits have x < 30, the TPC models in 
the PWN-determined ( strip are poor with x > 90. Po- 
larization swe ep data are limited fo r this pulsar. The pub- 
lished plots in I Johnston, et all (l2006|) give a « 34°, which is 
best matched with negative currents e < 0. These fits how- 
ever, imply a surprisingly large offset for the magnetic pole 
of 5(j) = —0.04 = 14°, which would put the magnetic axis at 
the leading edge of the radio pulse. Intriguingly, the best fit 
of all parameters is found with a OG PFF e = 0.05 model at 
C = 53°, a = 47° (x = 19). This very good light curve 
(Figure [9]) has the magnetic axis near the radio peak with 
5<j) = — 0.01, but would have a polarization sweep rate ~ 3x 
larger that that seen in present radio data. Again, improved 
radio polarization data might help shed light on the preferred 
inclination angle. 




The next pulsar, PSR B 1055-52 = PSR J1057-5226, is ap- 
preciably older. It shows a pulse that is nearly a square wave 
in the full energy band, alth ough the energy-dependent light 
curves of (lAbdo et all 12010) indicate that it is likely a tight 
double of separation A « 0.2 with strong bridge emission. 
The pulsar is too old to show a bright PWN, so no ( is avail- 
able from the X-rays. Luckily in the radio it has both a very 
broad pulse and interpulse and radio polarization fitting can 
strongly constrain both a and f. We adopt the values from 
(Welt evrede an d Wright 2009) and note the phase of the max- 
imum pol arization sweep is clearly at <j) = 0.08 (later) than 
plotted in (lAbdo et alll201Qh . Indeed, we find that the model 
fits prefer positive S(j), although the value is slightly larger 
than expected. The peaks on the model light curves are also 
stronger than in the data for these parameters; the match can 
be significantly improve d for models with fi nite 5w (not de- 
tailed here). As noted bv lWatters eTal (12009b . no viable TPC 
models lie anywhere near the RVM-determined angles. For 
OG there is a shallow fit minimum centered on the RVM an- 
gles. Superposed on this are x stripes, caused by the relatively 
coarse (25 bin) LAT light curve. With improved LAT data the 
fits should become more discriminating. 
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Fig. 8 . — Goodness of fit surface for PSR J1709-4429, panels as for Figure 
The X-ray determined £ is shown by the horizontal band. For this pulsar 
the smaller a RVM solution is clearly preferred. 



Fig. 9.— Light curves for two PFF models for PSR J1709-4429. The 
model at the RVM position has a rather weak second peak and requires a 
relatively large 5(f). A model at the PWN ( but slightly larger a provides the 
best fit at 86 ~ 0. 



For the TPC model there is a large region with modest x 
values at small a and with the best solution listed in Table 
1 . However these models are not only inconsistent with the ra- 
dio data, but all produce single pulses with one broad, approx- 
imately Gaussian, component and are unlikely to represent the 
true solution. If one allows any viewing angle, exceptionally 
good fits can in fact be found for OG with x slightly below 
the statistical minimum value. Despite an excellent match to 
the observed 7-ray pulse, this is also unlikely to represent the 
true solution. This underlines the fact that sensible solutions 
must be compatible with external angle constraints and, when 
these constraints are applied with care, one may use the multi- 
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wavelength data to rule out otherwise viable pulse models. 
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FIG. 10.— Goodness of fit surface for PSR J1057-5226, panels as for 
Figure|7] For this pulsar, interpulse emission allows a good RVM fit for both 
the inclination angle a and viewing angle £ (green circle). 

Finally, we discuss Geminga, the archetype of the old 7- 
selected pulsars, and the only one known prior to the launch 
of Fermi. This object does show an X-ray PWN, but like 
many old objects it is domina ted by bow shock structure 
(Pavl ov. Sanwal & Zavlin 2006) and so it will be very diffi- 
cult to extract torus/jet constraints on However we do know 
that despite very sensitive radio searches, no convincing de- 
tection of pulsed radio emission has been found, implying that 
|/3 1 = |C — ol\ is large. It has also been suggested that the 
thermal surface pulsations are most consistent with a nearly 
aligned rotator (small a, Pavlov, Sanwal & Zavlin 2006). 

Here, with no information from the radio phase we must 
allow the <fi to be freely fit. The results are interesting. For 
all B field structures, and in particular the PFF structure, the 
OG model only fits satisfactorily in a small region at small 
a, large (. This means that one is viewing a nearly aligned 
rotator from near the spin equator, i.e. very far from the mag- 
netic axis. This is in good agree ment with the expectations of 
(Romani and Yadigaroglul l 19951) and guarantees that the pul- 
sar should be radio faint. In figure Qj] we show surfaces com- 
puted for w = 0.18 and three trial e values. TPC shows best 
solutions over a range of angles and e. However, these have 
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Fig. 1 1 .— Goodness of fit surfaces for PSR J0633+1746 (Geminga), pan- 
els as for Figure [7] No quantitative external angle constraints are available, 
although the absence of radio emission suggests that (3 = ( — a > 15°. TPC 
shows a number of acceptable solutions; OG has few (albeit better) solutions 
at small a, large £. 
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Ax > 17 larger than the best OG solution and show only 
a single peak with bridge emission (Figure [12]). Thus while 
only few OG models are suitable, the light curve shape fits 
better and the viewing angles for Geminga are very well con- 
strained. It would be very interesting to measure £ and/or a 
from data at other wavelengths. 

In Table 1, we collect these model fit results, giving the 
values for the externally constrained fits for the unperturbed 
e = field geometry for both models. In addition we quote 
best fits close to the RVM angles and the best global fits, when 
appropriate. For Geminga, with no strong radio constraint we 
only list the global best fits for each model. 

3.4. Luminosities 

Now that we have estimated the corrections to the all- sky 
flux fa, it is worth checking how well the pulsars agree with 
the simple luminosity law Equation (1). A major challenge to 
this check is the distance imprecision; of the five pulsars con- 
sidered here only Vela and Geminga have parallax distance 
estimates. The other distances (derived from Dispersion Mea- 
sure DM) are subject to appreciable uncertainty. 

Geminga is in reasonable agreement with the rule only if 
it lies at the upper end of the present distance range. PSR 
J 105 7— 5226 is nominally somewhat more luminous than ex- 
pected, but can be accommodated if the DM distance is a 
slight overestimate. PSR J1709-4429, on the other hand, 
is much brighter than expected unless the true distance is 
~ 0.7 kpc or the flux correction factor is as small as fa ~ 
0.07 (or some combination of the two factors). Such changes 
would be quite surprising. It would be very interesting to 
obtain a parallax distance estimate for PSR J 1709— 4429 to 
eliminate that source of uncertainty. 

4. CONCLUSIONS AND FUTURE PROSPECTS 

In dWatters et alll2009|) we showed that patterns in the sep- 
aration of radio and 7-ray pulse components could be used 
to find acceptable geometrical parameters in dipole magne- 
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Fig. 12.— Light curves for two PFF models for PSR J0633+1746 
(Geminga). The left model is for the best OG solution, the right model for 
the best TPC solution away from f$ « 0. 



Fig. 13. — Pulsar 7-ray luminosity, corrected from the observed phase- 
averaged flux to all-sky vs. the spindown luminosity. OG/PFF/e = points. 

tosphere models for pulsar emission. Here we show that 
with high quality LAT light curves and reference to multi- 
wavelength information on the inclination and viewing angles 
the direct comparison between the model light curve shape is 
even more constraining. Within the restricted context of the 
present computations - dipole fields, simple expressions for 
the gap ranges and 7-ray emissivity on the radiating surface - 
we can already make some strong statements. 

• Emission starting above the null charge surface (OG model) 
is strongly statistically preferred over models which have sub- 
stantial emission starting from the stellar surface (eg. the TPC 
model). For PSR J1952+3252 the fit for the TPC model is not 
much worse than that for OG, but for none of the pulsars mod- 
eled here is it better, if external angle constraints are used. 

• More realistic dipole field geometries, including sweep back 
and fields enforcing co-rotation produce improved matches to 
the pulsar light curves when one focuses on solutions near 
the orientation angles a and £ implied by lower energy data. 
Some additional improvement can be found by using test 
models for field-induced current perturbations, but these are 
evidently not yet sufficiently realistic for detailed fits. 

• The correction from observed flux to pulsar luminosity 
can be substantial. With these corrections a simple constant 
<£, L 1 oc E(oc current) law provides an improved descrip- 
tion of the data. Nevertheless, the large departure for PSR 
J 1709— 4429 suggests additional factors are needed to under- 
stand pulsar luminosities. 

• The best-fit global models are often several degrees off of 
the externally known ( and a. This difference is statistically 
significant and implies that the LAT data have the power to 
reveal perturbations to the simplified field geometry. 

The robust success of the pulse profile computations de- 
scribed here for angles near those available from other data 
make it likely that the true magneto spheric geometry is fairly 
close to our simplified PFF model. It must be re-emphasized 
that this PFF model is not a complete physical description - 
and indeed the whole electrodynamic picture of 'Holloway- 
type' outer gaps is still of questionable validity. How- 
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TABLE 1 

PFF Model Angles and Efficiencies 



Name 


\og(E) 


w 


d a 


a/C a/C 


a/C/e 




fn 


a/C/e 




fn ^>0.1GeV a VOG 




ergs -1 






i — External — > 




TPC 






OG — 





Vela 


36.84 


0.012 


J1952+3252 


36.57 


0.016 


J1709-4429 


36.53 


0.017 


Geminga 


34.52 


0.175 


J1057-5226 


34.48 


0.183 



0-287±S:gi? 

2.0±0.5 
1.4-3.6 

u - zo -0.06 

0.71±0.2 



best 
best 

b 

best 



in reg 
global 

b 



71/84/0 
86/66/-0.2 



in reg 

34/53 108/127[72/53] 31/49/0 
best global 

50/26/0 

105/1 11 [75/69] 75/69/0 
best global 41/08/-0.2 



201/+0.02 1.06 72/64/0 

71/64/+0.2 
82/71/-0.2 
19/+0.05 1.10 66/78/0 
18/+0.05 1.10 71/71/-0.2 
147/+0.10 1.30 36/56/0 

47/53/+0.2 
258/+0.41 0.93 21/89/0 

414/+0.12 0.91 75/69/0 
37/+0.03 0.38 7/87/-0.2 



120/-0.03 1.03 

88/-0.02 0.98 

34/-0.03 0.86 

14/+0.05 0.84 

14/+0.05 0.70 

30/-0.04 0.89 

19/-0.01 0.76 

241/+0.87 0.13 

106/+0.11 0.82 

5/+0.11 0.09 



879.4±5.4 

13.4±0.9 
124±2.6 
338.1±3.5 
27.2±1.0 



0.015 
0.24 
0.05 
0.45 



a From Abdo et al ( 2010), flux units 10 -11 ergcm 2 s -1 
b d^/d(j)\max « +3.5°/deg-Weisbergetal. (1999). 

ever, other scenarios need to match the light curves quite 
closely to be competitive - this likely forces them into rather 
similar magnetosphere geometries. Comparison with other 
approaches, such as numerical magnetosphere simulations, 
should be fruitful. With more pulsar light curves and higher 
quality data being provided by the Fermi LAT, we should be 
able to use the departures from the predictions of these simple 
beaming laws to infer geometrical modifications (and work 
toward their physical origin) and to follow the evolution of 
the gap geometry and radiation with pulsar age. In a follow- 
ing paper we describe the additional constraints on w(E), Sw, 



the radial extent of the emission zone, and the pulsar evolu- 
tion with age that can be extracted from the full LAT pulsar 
data set. This should set the context for true radiation mod- 
els, which can take the next step in revealing the underlying 
magnetospheric physics. 
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APPENDIX 

For the convenience of other researchers we present here figures that summarize the pulse properties of the TPC and OG models 
for the PFF dipole geometry assumptions. These may be compared with the corresponding figures in Watters et al. (2009). 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



a=15° a = 35° a = 55° a = 75° a = 90° 

Fig. 14. — Light curves for the Two Pole Caustic (TPC) model. Each panel shows curves for four values of the gap width w. The curves are labeled with the 
number of major peaks and the peak separation, in percent. 
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20 40 60 80 20 40 60 80 20 40 60 80 



Magnetic Inclination a 

Fig. 16. — Pulse properties in the a — £ plane for two pulsar emission geometries using the PFF field. Contours show the separation A of the principal 7-ray 
pulse peaks. Bold lines mark A at 0.1 intervals (key bottom middle panel); fine lines mark 0.02 intervals. The background shading (key bottom left panel) gives 
the flux correction factor /q; large values mean that for the given a and the pulsar displays less flux than the sky average. White regions mean little or no 
flux from the modeled gaps. The dashed diagonal band in each panel indicates the region where lower altitude radio pulsations are likely detected. Objects seen 
outside these bands tend to be Geminga-like radio-quiet pulsars. 



